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Thermal radiation from hot polar cap in pulsars 
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Abstract. Thermal radiation from hot polar caps is ex- 
amined in radio pulsars with drifting subpulses. It is ar- 
gued that if these subpulses correspond to sparking dis- 
charges of the inner acceleration region right above the 
polar cap surface then a simple relationship between the 
observed subpulse drift rate in radio and thermal X-ray 
luminosity from the polar cap heated by sparks should 
exist. This relationship is derived and tested in pulsars 
for which an appropriate good quality data is available. 



1. Introduction 

Almost 40 years have passed since the discovery of pul- 
sars and the mechanism of their coherent radio emission is 
still not known. The puzzling phenomenon of drifting sub- 
pulses is widely regarded as a powerful tool for the investi- 
gation of mechanisms of pulsar radio emission. In the clas- 
sical model of Rudcrman & Sutherland (1975; RS75 hence- 
forth) the spark-associated subbeams of subpulse emission 
circulate around the magnetic axis due to E x B drift 
of spark plasma filaments. This model is widely regarded 
as the most natural and plausible explanation of drifting 
subpulse phenomenon, at least qualitatively. Despite its 
popularity the RS75 model is known to suffer from the 
so-called binding energy problem (for review see Usov &i 
Melrose 1995, 1996). In fact, the cohesive energy of surface 
iron ions were largely overestimated in RS75 and the "vac- 
uum gap" envisioned by RS75 was impossible to form. A 
number of attempts have been made to resolve this prob- 
lem, including a quite exotic proposal that radio pulsars 
with drifting subpulses are bare polar cap strange stars 
(BPCSS) rather than neutron stars (Xu, Qiao & Zhang, 
1999). Gil & Melikidze, (2002) argued that the formation 
of RS75 gap above the neutron star polar cap was, in prin- 
ciple, possible, although it required a very strong surface 
magnetic field, much stronger than the dipolar component 
inferred from the observed spin-down rate. Growing evi- 
dence of such strong non-dipolar surface field accumulates 
in the literature, both observational and theoretical (see 
Urpin & Gil 2004, for short review). 



Even if the RS75 gap was possible to form, it would 
not automatically solve the mystery of drifting subpulses. 
It is well known that the original model of RS75 predicts 
too fast a drifting rate (e.g. Deshpandc & Rankin, 1999; 
DR99 henceforth). Motivated by this issue Gil, Melikidze 
& Geppert (2003; GMG03 henceforth) developed further 
the idea of the inner acceleration region above the polar 
cap by including the partial screening due to thermionic 
ions flow from the surface heated by sparks. We will call 
this kind of the inner acceleration region the "partially 
screened gap" (PSG henceforth). The slow drift corre- 
sponds to E x B drift of spark generated electron-positron 
pairs, until the total charge density within the partially 
screened gap reaches the co-rotational value. Since the 
PSG potential drop is much lower than in the RS75 model, 
the intrinsic drift rate is compatible with observations (for 
details sec GMG03). 

A distinguishing property of PSG model is relatively 
high predicted heating rate of the polar cap surface, com- 
patible with observations (in contrast to RS75 gap, which 
overheated the polar cap). The space charge limited model 
(Arons & Scharlcman, 1979, AS79 henceforth; Zhang & 
Harding, 2000; Harding & Muslimov, 2002, HM02 hence- 
forth) predicts a much lower polar cap heating rate. On 
the other hand, in the BPCSS model (Xu, Qiao, & Zhang, 
1999) no hot spot is expected due to the high thermal con- 
ductivity at the BSS surface. Thus, measuring the ther- 
mal X-ray luminosity from heated polar caps can poten- 
tially reveal the nature of the inner acceleration region in 
pulsars. This can help us to understand a mechanism of 
drifting subpulses, which appears to be a common phe- 
nomenon in radio pulsars. Recently Wcltevrede, Edwards, 
& Stappcrs (2006, WES06 henceforth) presented results 
of a systematic, unbiased search and found that the frac- 
tion of pulsars showing drifting subpulses is at least 55 %. 
They concluded that the conditions for drifting mecha- 
nism to work cannot be very different from the emission 
mechanism of radio pulsars. 

In order to test different available models of inner ac- 
celeration region in pulsars Zhang, Sanwal & Pavlov 2005; 
ZSP05 henceforth) observed the best studied drifting sub- 
pulse radio pulsar PSR B0943+10 with the XMM-Newton 
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observatory. Their observations were consistent with PSG 
formed in strong, non-dipolar magnetic field just above 
the surface of very small polar cap. Recently, Kargaltsev, 
Pavlov & Garmire (2006, KPG06 henceforth) observed the 
X-ray emission from the nearby PSR B1133+16 and found 
that this case is also consistent with the thermal radiation 
from a small hot spot (again much smaller than the canon- 
ical polar cap). PSR B1133+16 is almost a twin of PSR 
B0943+10 in terms of P and P values and, interestingly, 
both these pulsars have very similar X-ray signatures, in 
agreement with our PSG model (see Table 1). 

In this paper we generalized the treatment of ZSP05 
and developed detailed model for thermal X-ray emission 
from radio drifting pulsars. The model matches the ob- 
servations of PSRs B0943+10 and B1133+16 well, both 
in radio and X-rays. There is a number of pulsars with 
measured thermal X-ray radiation from small hot polar 
cap but unfortunately their drifting properties are not yet 
known. This is likely to change in the near future due to 
new sophisticated methods for analysis of intensity fluctu- 
ations in weak pulsars being developed (e.g. WRS06 and 
references therein). 

2. PSG model of the inner acceleration region 

As already mentioned, growing evidence appears, both 
observationally and theoretically, that the actual surface 
magnetic field B s is highly non-dipolar. Its magnitude 
can be described in the form B s = bBd (Gil & Sendyk 
2000; GS00 henceforth), where the enhancement coeffi- 
cient b > 1 and B d = 2 x 10 12 (PP_i 5 ) 1/2 G is the canon- 
ical, star centered dipolar magnetic field, P is the pulsar 
period and P_is = P/10" 15 is the period derivative. 

The polar cap is defined as the locus of magnetic field 
lines that penetrate the so-called light cylinder. Conven- 
tionally, the polar cap radius r pc = 1.45 x 10 4 P~ 05 cm, 
and its surface area A pc ~ 2 x 10 8 P _1 cm 2 . In the 
case of non-dipolar surface field the polar cap area must 
shrink due to flux conservation of the open field lines 
A pc Bd = A p B s . Thus, the surface area of the actual polar 
cap A p = b A pc , regardless of the actual shape of the 
polar cap. Consequently, one can write the polar cap ra- 
dius in the form r p = 6~°' 5 r pc (GS00), realizing however 
that this is only a characteristic dimension of the actual 
polar cap. First, the canonical expression for the radius 
of the polar cap (GJ69) is only a geometrical approxima- 
tion (e.g. Michel 1973). Secondly, the presence of strong 
surface magnetic field anomalies should make the shape 
of the polar cap quite irregular. 

The charge depleted inner acceleration region above 
the polar cap results from the deviation of a local charge 
density p from the co-rotational charge density (Goldrc- 
ich & Julian 1969) p GJ = —ft ■ B s /2ttc w B s /cP. For 
isolated neutron stars one might expect the surface to con- 
sist mainly of iron formed at the neutron star's birth (e.g. 
Lai 2001). Therefore, the charge depletion above the po- 



lar cap can result from bounding of the positive IJjFe ions 
(at least partially) in the neutron star surface. As demon- 
strated by recent exact calculations of Mcdin & Lai (2006) , 
the iron chains are strongly bound in magnetic field close 
to 10 14 G. This can only occur if the surface magnetic 
field in actual pulsars is dominated by strong non-dipolar 
components. If this is really possible, then due to signif- 
icant bounding the positive charges cannot be supplied 
at the rate that would compensate the inertial outflow 
through the light cylinder. As a result, a significant part 
of the unipolar potential drop develops above the polar 
cap, which can accelerate charged particles to relativistic 
energies and power the pulsar radiation mechanism. 

The ignition of cascading production of electron- 
positron plasma is crucial for limitation of growing gap 
potential drop. The accelerated positrons will leave the 
acceleration region, while the electrons will bombard the 
polar cap surface, causing a thermal ejection of ions. This 
thermal ejection will cause partial screening of the accel- 
eration potential drop AV corresponding to a shielding 
factor 77 = 1 — Pi/pGJ (see GMG03 for details), where pi 
is charge density of ejected ions, AV = / cP)B s h 2 
is the potential drop and h is the height of the accelera- 
tion region. The gap potential drop is completely screened 
when the total charge density p = pi -\- p+ reaches the co- 
rotational value Pgj- 

GMG03 argued that the actual potential drop AV 
should be thermostatically regulated and the quasi- 
equilibrium state should be established, in which heating 
due to electron bombardment is balanced by cooling due 
to thermal radiation. The quasi-equilibrium condition is 
Qcooi = Qheat, where Q coo i = crT^ is a cooling power 
surface density by thermal radiation from the polar cap 
surface and Qheat = ~/m e c 3 n is heating power surface den- 
sity due to back-flow bombardment, 7 = eAV/m e c is 
the Lorentz factor, n = uqj — n\ — tjugj is the number 
density of back-flowing plasma particles depositing their 
kinetic energy at the polar cap surface, r\ is the shielding 
factor, m is the charge number density of thermionic ions 
and n GJ = p GJ /e = 1.4 x 10 11 6P° 1 5 5 P-°- 5 cm- 3 is the 
corotational charge number density. It is straightforward 
to obtain an expression for the quasi-equilibrium surface 
temperature in the form 

r s = (6.2 x io^xp^s/p) 1 / 4 ?? 172 ^ 172 - (1) 

Let us now interrelate the accelerating potential drop 
AV and the perpendicular (with respect of the magnetic 
field lines) electric field AP which causes E x B drift. Fol- 
lowing the original method of RS75 we can argue that 
the tangent electric field is strong only at the polar cap 
boundary where AP = 0.5AV/h — ^(n / 'cP)B s h (sec Ap- 
pendix A in GMG03 for details). Due to the E x B drift 
the discharge plasma performs a slow circumferential mo- 
tion with velocity Vd = cAE/B s = ijTrh/P. The time in- 
terval to make one full revolution around the polar cap 
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boundary is P3 s» 2irr p /vd- One then has 
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If the plasma above the polar cap is fragmented into fil- 
aments (sparks) which determine the intensity structure 
of the instantaneous pulsar radio beam, then in princi- 
ple, the tertiary periodicity P3 can be measured/estimated 
from the pattern of the observed drifting subpulscs (Desh- 
pande & Rankin 1999, Gil & Sendyk 2003). According to 
RS75, P3 = NP 3 , where N is the number of sparks con- 
tributing to the drifting subpulse phenomenon observed 
in a given pulsar and P3 is the primary drift periodicity 
(distance between the observed subpulse drift bands). On 
the other hand N ss 2-nr v j2h (GS00). Thus, one can write 
the shielding factor in the form 77 pb (l/27r)(P/P 3 ) , which 
depends only on an easy-to-measure primary drift peri- 
odicity. Apparently, the shielding parameter 77 should be 
much smaller than unity. 

The X-ray thermal luminosity is L x 
1.2 x l0 32 (P- 15 /P 3 )(rih/r p ) 2 erg/s, which can be com- 
pared with the spin-down power E = IflCl = 3.95/45 x 
10 31 P_i 5 /P 3 Cr g/ S; where I = Ii 5 10 45 g cm 2 is the neu- 
tron star moment of inertia (bellow we assume that 
Z45 = 1). Using equation ([2]) we can derive the ther- 
mal X-ray luminosity and its efficiency as L x — 2.5 x 
10 31 (P_i5/P 3 )(P/p3) 2 , or m the simpler form represent- 
ing the efficiency with respect to the spin-down power 

f=°- 63 (£> 2 - < 3 > 

which is very useful for comparison with observations. One 
should realize that this equation holds only for thermal 
X-rays from hot spot and cannot be applied neither to 
cooler radiation from the entire stellar surface nor to the 
magnctosphcric component. 

We can see that L x in these equations depends only 
on radio observables. It is particularly interesting and im- 
portant that both equations above do not depend on de- 
tails of the sparking gap model (r],b,h). Although one 
has to be careful whether all our assumptions are satis- 
fied in real pulsars, it seems that we have found a very 
useful, relatively easy testable (at least for the order-of- 
magnitude) relationship between the properties of drift- 
ing subpulses observed in radio band and the charac- 
teristics of thermal X-ray emission from the polar cap 
heated by sparks associated with these subpulscs. For 
PSRs B0943+10 and B1133+16, which are the only two 
pulsars in which both P 3 and L x are measured/estimated 
(Table 1), the above equation holds very well. Interest- 
ingly, in both cases L x /E ~ 10~ 3 . 

Using equation ([2]) we can write the polar cap temper- 
ature in the form 

-1/2 

T s = (5.1 x 10 6 K)6 1 / 4 py i 4 5 P" 1 / 2 [ ^ ] , (4) 




A p is the actual emitting sur- 
face area (bolometric) . Since A bo i can be determined from 
the black-body fit to the spectrum of the observed hot- 
spot thermal X-ray emission, the above equations can be 
regarded as independent of details of the sparking gap 
model and depending only on combined radio and X-ray 
data, similarly as in equations Q. 

3. Comparison with observational data 

Table 1 presents the observational data and predicted val- 
ues of a number of quantities for two pulsars, which we be- 
lieve show clear evidence of thermal X-ray emission from 
the spark heated polar caps as well as they have known 
values of tertiary subpulse drift periodicity. The predicted 
value of P3 and/or L x were computed from cq.(4), while 
the predicted values of T s were computed from cq.(5), with 
b = Ap C /Aboi determined observationally. The actual sur- 
face magnetic field B s = bB<i, where Bd is the canonical 
dipolc magnetic field at the polar cap. As a principle, the 
predicted (Prcd) values were obtained from observational 
(Obs) values through Eq.(3). 

PSR B0943+10. This is the best studied drifting 
subpulse radio pulsars with P = 1.09 s, P-15 = 3.52, 
E = 10 32 erg s"\ P 3 = 1.86P, P 3 = 37.4P and N = 
P3/P3 = 20 (DR99). It was observed by ZSP05, who ob- 
tained an acceptable thermal BB fit with bolometric lu- 
minosity L x = (5±?;|) x 10 28 erg s _1 and thus L x /E = 
(0-49±g'?|) x 10" 3 . The bolometric polar cap surface area 
A bol = 10 7 [T s /(3 x 10 6 K)r 4 cm 2 ~ (1± 4 °) x 10 7 cm 2 
is much smaller than the conventional polar cap area 
A pc = 6 x 10 s cm 2 . This all correspond to the best fit tem- 
perature T s ~ 3.1 x 10 6 K (see Fig. 1 in ZSP05). The pre- 
dicted value of L x /E calculated from equation agrees 
very well with the observational data. The surface temper- 
ature T s calculated from equation (J4j) with b = A pc /Aboi 
is also in good agreement with the best fit. The shield- 
ing factor rj — 0.09, thus more than 90 % of the available 
vacuum potential drop is screened. 

PSR B1133+16. This pulsar with P = 1.19 s, 
P_i5 = 3.7, and E = 9 x 10 31 erg s _1 is almost a twin of 
PSR B0943+10. KPG06 observed this pulsar with Chan 
dra and found an acceptable BB fit L x /E = (0.77^ 



+0.13N 
15) 



10~ 3 , A bo i = (0-5±g;|) x 10 7 cm 2 and T s w 2.8 x 
10 6 K. These values are also very close to those of PSR 
B0943+10, as should be expected for twins. Using equa- 
tion © we can predict P 3 /P = 27±2 for B1133+16. In- 
terestingly, Nowakowski (1996) obtained fluctuation spec- 
trum for this pulsar with clearly detected long period fea- 
ture corresponding to about 32P. Most recently, WES06 
found P3 / P = 3 ± 2 and long period feature corresponding 
to (33±3)P in the fluctuation spectrum of PSR Bl 133+16. 
The latter value seem to coincide with that of Nowakowski 
(1996), as well as with our predicted range of P3. We there- 
fore claim that this is the actual tertiary periodicity in 
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Name 


P3/P 


Ps/P 


LJE 


x 10 3 


b 


rp(obs) 


y(prcd) 


B d 


B s 


PSR B 


Obs. 


Obs. Pred. 


Obs. 


Pred. 


Ape/ ^bol 


10 b K 


10 b K 


10 i2 G 


10 i4 G 


0943 + 10 
1133+16 


1.86 

q+2 

J -2 


37.4 36+^ 
(33±1) 27tl 


49+^6 

77+ - 13 
u -' '-0.15 


0.45 


60+ 140 

ou -48 
11 1+16-6 
L1 - L -5.6 


q 1 +0.9 
°- ± -l.l 

2.811:1 


3.3±i-« 
9 1+0.5 


3.95 
4.25 


9 07+5.53 

z -°'-1.90 

47+ 71 



PSR B1133+16 and show it in paranthcsis in Table 1. 
It is worth noting that P3/P3 = 33 ± 3 is quite close to 
37.4 measured in the radio twin PSR B0943+10. Note also 
that the number of sparks predicted from our hypothesis 
is N = P3/P3 = (33 ± 3)/(3 ± 2) = lltf, so it can also 
be close to 20, as in the case of twin PSR B0943+10. The 
shielding factor r\ = 0.051q'q2. Thus, only several % of the 
vacuum gap potential drop is available for acceleration. 

4. Conclusions and discussion 

Within the partially screened gap model of the inner ac- 
celeration region in pulsars developed by GMG03 we de- 
rived a simple relationship between the X-ray luminosity 
L x from the polar cap heated by sparks and the tertiary 
periodicity P3 of the spark-associated subpulse drift ob- 
served in radio band. In PSRs B0943+10 and B1133+16 
for which both L x and P3 arc known, the predicted rela- 
tionship between observational quantities holds quite well. 

Both the heating and the drifting rate depends on 
the gap potential drop. In this paper we made the point 
that there is continuum of cases between pure vacuum 
gap and the space charge limited flow. The original RS75 
model predicts much too high a subpulse drift rate and 
an X-ray luminosity. Other available acceleration mod- 
els predict too low a luminosity and the explanation of 
drifting subpulse phenomenon is generally not clear (see 
ZSP05 for more detailed discussion). Approximately, the 
bolometric X-ray luminosity for the space charge limited 
flow (Arons & Sharlcman 1979) is about (10~ 4 + 10~ b )E 
(Harding & Muslimov 2002), and for the pure vacuum 
gap (RS75) is about (10" 1 +- 10" 2 )P (ZSP05), while for 
the partially screened gap (GMG03) is ~ 10~ 3 _E (this pa- 
per). The latter model also predicts right E x B plasma 
drift rate. Thus, combined radio and X-ray data are con- 
sistent only with the partially screened gap model, which 
requires very strong (generally non-dipolar) surface mag- 
netic fields. Observations of the hot-spot thermal radiation 
almost always indicate bolometric polar cap radius much 
smaller than the canonical value (by an order of magni- 
tude). Most probably such a significant reduction of the 
polar cap size is caused by the flux conservation of the 
non-dipolar surface magnetic fields connecting with the 
open dipolar magnetic field lines at distances much larger 
than the neutron star radius. A small part of this reduc- 
tion can follow from the fact that discrete sparks do not 
heat the entire polar cap area (say about a half), but un- 



certainties in determination the size of the polar cap are 
likely to account for this effect. 
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